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ABSTRACT 

We present a new population of z > 2 dust-reddened, Type 1 quasars with 0.5 < 
E{B — V) < 1.5, selected using near infra-red (NIR) imaging data from the UKIDSS- 
LAS, ESO-VHS and WISE surveys. NIR spectra obtained using the Very Large Telescope 
(VLT) for 24 new objects bring our total sample of spectroscopically confirmed hyperlu- 
minous (>10^^Lq), high-redshift dusty quasars to 38. There is no evidence for reddened 
quasars having significantly different Ha equivalent widths relative to unobscured quasars. 
The average black-hole masses (^1O®-1O^*’M0) and bolometric luminosities (^10^^ erg/s) 
are comparable to the most luminous unobscured quasars at the same redshift, but with a 
tail extending to very high luminosities of erg/s. Sixty-six per cent of the reddened 

quasars are detected at > 3cr at 22fim by WISE. The average 6fim rest-frame luminosity is 
logio(L6;im/erg s”^)=47.1±0.4, making the objects among the mid-infrared brightest AGN 
currently known. The extinction-corrected space-density estimate now extends over three 
magnitudes (—30 < Mi < —27) and demonstrates that the reddened quasar luminosity func¬ 
tion is significantly flatter than that of the unobscured quasar population at z = 2 — 3. At 
the brightest magnitudes. Mi < —29, the space density of our dust-reddened population 
exceeds that of unobscured quasars. A model where the probability that a quasar becomes 
dust-reddened increases at high luminosity is consistent with the observations and such a de¬ 
pendence could be explained by an increase in luminosity and extinction during AGN-fuelling 
phases. The properties of our obscured Type 1 quasars are distinct from the heavily obscured, 
Compton-thick AGN that have been identified at much fainter luminosities and we conclude 
that they likely correspond to a brief evolutionary phase in massive galaxy formation. 

Key words: galaxies:active, (galaxies:) quasars: emission lines, (galaxies:) quasars: general, 
(galaxies:) quasars: individual 


1 INTRODUCTION 

Luminous dust-obscured quasars have been postulated to repre¬ 
sent the missing evolutionary link between merger -induced star- 
burst galaxies and ultraviolet (UV) bright quasars jSanders_e£^ 
198^ in theoretical models of galaxy formation (e.g. lHonkins et al' 


20081: iNaravanan 

et alJl201Gh. Samnles of obscured active galac- 

tic nuclei (AGN) have been selected both in the hard X-ray (e.g. 

Brusa et al.ll201(]l: 

Civano et alj2012l: LaMassa et al.ll2013h and the 

mid-infrared (e.g. 

Stem et al.l2012l: Assef et al.l2013l; Donlev et all 
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[ 20 II . However, these surveys do not distinguish between a cen¬ 
trally concentrated obscuring medium, on parsec scales, close to 
the b lack-hole accretion disk dUrrv & Padovanilll995l : lAntonuccil 
Il993h . and larger-scale obscuration due to gas and du st clouds in 
the quasar host galaxy jMartmez-Sansigre et alJl^OSl) . Most hard 
X-ray and mid-infrared surveys (e.g. with the Spilzer Space Tele¬ 
scope) have covered relatively small areas of sky, probing obscured 
black-hole growth in relatively low-luminosity AGN, although the 
situation is now changing at m id infrared wavelen gths with data 
from the WISE All-Sky Survey iWright et alj|201(]h . In these faint 
samples, the number of obscured AGN is found to decrease with 
increasing luminosity. Constraints on obscured black-hole growth 
at the highest quasar luminosities however remain highly uncertain 
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and the relationship between hyperluminous (> lO^^L©) obscured 
quasars and the now well-studied obscured AGN at low luminosi¬ 
ties is also uncertain. 

X-ray surveys also show a curious lack of high Eddington 
ratio (L/LEdd > 0-1) objects with Nh between 3x10^^ and 
3x10^^ cm“^ ( lFabianll2012l : lRaimundo et alj|2010l) . Objects with 
such properties, i.e. high luminosity combined with a significant 
gas/dust column, could possess massive AGN-driven outflows. The 
apparent upper limit to the distribution of Eddington ratios could be 
a selection bias, as most X-ray surveys do not cover a large enough 
area to detect the rarest and brightest quasars with high Edding¬ 
ton ratios. Such objects have traditionally been selected through 
their rest-frame UV emission which is very sensitive to dust extinc¬ 
tion and therefore biased against detecting luminous quasars with 
even moderate gas column densities. The Sloan Digital Sky Sur¬ 
vey ( SPSS) Data Release (DR) 7 quasar sample jSchneider et al.l 
120101) now mitigates the situation to a degree. The flux limit in 
the redder t-band allows some intrinsicall y very luminous ob - 


l y V 

jects with E{B — V) ~0.5 to be included ^Richards et al.l| 2003h. 


The new SDSS-III BOSS DRIO quasar survey Jparis et ak boi^i) 
reaches fainter optical magnitudes than the SPSS and has signifi¬ 
cantly increased the number of quasars known at the main epoch 
of galaxy formation at z ~ 2 — 3. Populatio ns of red quasars 
are detected in the BOSS data ( IRoss et al.l20l4) including the first 
large sample of narrow-line T ype 2 obscured quasars at high red- 
shifts dAlexandroff et al.ll2013h . so me of which also sho w evidence 
for broad-lines in the near infrared dOreene et alj|201'i) . However, 
even with the fainter flux-limit, BOSS cannot find quasars with 
E{B -V) >0.6. 


Wide-field infrared surveys now cover much of the observable 
sky, providing the sensitivity to detect dusty, intrinsically luminous, 
broad-line quasars. The 2MASS survey has i dentified optically ob¬ 
scured quasars that are also radio-loud (e.g. iGlikman et alJlrOO^ : 
lUrrutia et alJl2009l : lGlikman et al.ll2012h : the radio selection using 
FIRST data is necessary to eliminate contaminant populations but 
means that the majority of AGN are excluded by the selection. In 
radio-loud red quasars (e.g. IWebster et al.lll995h . the red colours 
can a rise as a result of the presence of a red synchrotron compo¬ 
nent dFrancis et al.ll200ll : IWhitine et alj|200i]) rather than due to 
dust. The flux-limit of the 2MASS survey has also meant that these 
red quasars are predominantly at low-redshift altho ugh with some 
of the most luminous quasars seen out to z > 2 dGlikman et al.l 
[201^ . A concern at low-redshifts is that host galaxy contamina¬ 
tion can redden the observed colours of these quasars. The galaxy 
spectral energy distribution (SED) peaks at ~l/im, corresponding 
to the inflection point in the AGN SED. This feature is in the NIR 
JHK filters at z ~ 0.2 — 1. However, at z -^0.7, the FIRST- 
2MA SS red quasars are s een to be associated with merging galaxy 
hosts dUrrutia et al.ll200^ . consistent with an interpretation where 
the objects are in a transition phase from starburst to unobscured 
quasar. Extending these investigations to focus on the main epoch 
of galaxy and black-hole formation at z ~ 2 — 3 will be in in¬ 
structive in understanding the importance of this obscured phase in 
galaxy formation. 

The WISE All-Sky Survey operating at mid-infrared wave¬ 
lengths, has recently led to the identification of a new population of 
hyperluminous dust-obscured galaxies (HyLIRGs), whose spectral 
energy distribution s (SEDs) are consistent with a dominan t contri¬ 
bution from AGN dEisenhardt et alj|2012l : IWu et alj|2012h . These 
HyLIRGs are also postulated to represent the evolutionary link be¬ 
tween massive starbursts and luminous quasars. The WISE samples 
are complementary to near infrared selected samples as they enable 


the selection of much more highly obscured AGN. The AGN that 
have so far been identified by WISE are seen to be Compton thick 
dStem et al.l2014h . Broad emission lines originating from gas close 
to the black-hole are therefore completely concealed and estimating 
virial black-hole masses is no longer possible. 

The wide-fi eld N IR UKIDSS Large Area Survey (LAS) 
2007h and VISTA Hemisphere Survey (VHS; 
allow us to probe much fainter than was pos¬ 
sible using 2MASS and enable selection of complementary sam¬ 
ples of obscured quasars to the WISE survey. In a survey over 
> 1000 deg^ of the UKIDSS and VHS data we were a Me to identify 
a population of z 2 — 3, E{B — V)^ 0.8 quasars dBanerii et al.l 
l2012Ll20l^ - hereafter B12 and B13 respectively), demonstrating 
that contaminant populations can be eliminated without having to 
rely on radio pre-selection. Our new quasars are among the most 
bolometrically luminous quasars known (Lboi >10'*^ erg/s) but, 
as a result of their significant dust-extinctions, are often invisible 
in wide-field optical surveys like SDSS. The sample represents a 
population of hyperluminous quasars at z 2 — 3 whose extinc¬ 
tion properties are intermediate between UV-selected, essentially 
unobscured, quasars and the very heavily obscured AGN emerging 
from WISE. Similar obscured luminous quasars showing strong ev¬ 
idence for corresponding to the AGN feedback phase, are also now 
emerging fro m multi-wavelength degree-scale galaxy surveys such 
as COSMOS jBrusa et al.l201()l : lBongiomo et alj2014l : lBrusa et al.l 
l2014h . Indeed the brightest obscured quasar selected from the COS¬ 
MOS survey is the faintest obscured quasar selected in our wide- 
field UKIDSS search for similar objects (B12). The ability to iso¬ 
late the luminous AGN population with intermediate levels of ex¬ 
tinction means that it is now possible to investigate objects with 
extinctions corresponding to reddenings oi 0 < E{B — V) < i. 
New insight into the importance of non-spherical geometries (i.e. 
‘unified schemes’) and evolution, involving specific phases during 
AGN lifetimes, for explaining the distribution of extinction proper¬ 
ties can be expected. 

Our survey also led to the identification of the prototype 
quasar; ULASJ1234-1-0907, the reddest broad-line quasar currently 
known. With an Eddington ratio of 0.65 (B12) and a neutra l 
hydrogen column-density of 9x10^^ cm~^ dBanerii et al.ll20l3) . 
ULASJ 1234-1-0907 lies in the sparsley-populated high L /LBdd and 
intermediate gas column regime found from the statistics of hard 
X-ray sources. Among the most luminous quasars known at both 
X-ray and far-infrared wavelengths, ULASJl234-1-0907 may repre¬ 
sent the em erging phase of a su permassive black-hole shortly after 
a starburst dBanerii et al.ll2014h . Larger surveys for NIR-selected, 
dust-obscured, luminous quasars may unearth more objects like 
ULASJl 234-1-0907. 

The current paper presents a much larger survey for heavily 
dust-reddened, hyperluminous broad-line quasars at z = 2 — 3, 
where the abundance of gas-rich, star-forming galaxies provides 
a natural explanation for the presence of extended obscuring me¬ 
dia in these luminous quasars. Combined with objects presented in 
B12 and B13, the new sample allows us to characterise in detail 
the demographics of the luminous reddened quasar population at 
z 2.5, including the direct calculation of the space density over 
a three-magnitude interval. We emphasise that the parameter space 
sampled by our survey is very different from that of X-ray surveys 
for obscured AGN in deep fields covering only a few square de¬ 
grees. A companion paper (Alaghband-Zadeh et ah, in preparation) 
investigates the spatially resolved properties of the new reddened 
quasars from the IFU data presented here. 

The structure of the paper is as follows. First, the observational 
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data, including the photometric selection (Section 2) and spectro¬ 
scopic follow-up (Section 3), is described. The generation of the 
key physical properties of the quasars, including extinction esti¬ 
mates, luminosities and black-hole masses, is covered in Section 4, 
prior to consideration of the object properties derived from WISE 
photometric data in Section 5. The primary scientific result, the cal¬ 
culation of the population space density and luminosity function is 
presented in Section 6, followed by a summary of the main results 
in Section 7. Throughout the paper we assume a flat ACDM cos¬ 
mology with ho=0.70, fim = 0.30 and Ha= 0.70. The VHS mag¬ 
nitudes are on the AB system where the conversions from Vega to 
AB used are: J = -bO.937, H = -bl.384, K = 1.839. All WISE 
magnitudes and colours are on the Vega system. 

2 PHOTOMETRIC SELECTION 

2.1 Colour Selection of Reddened Quasars 

We make use of new wide-held infrared photometry over ~1500 
deg^ from VHS and WISE, in order to identify luminous reddened 
quasar candidates in the southern hemisphere. In addition, we also 
use >3000 deg^ of imaging data fr om the UKIDSS Large Area Sur¬ 
vey (LAS) ( [Lawrence et al.ll2007h in the northern hemisphere, to 
select a sample of t he reddest quasars, such as ULASJ1234-1-0907 
dBanerii et al.ll201^ . The details of the colour selection has been 
presented in B13 but we provide a summary here. Candidates must 
be morphologically classihed as point sources in the Tf-band (i.e. 
with kclass = — l: lGonzalez-Solares et al]|201 ih and possess ex¬ 
tremely red NIR colours of (J—TfjAB > 1.6. The morphology cut 
is necessary to ensure that we isolate a population of quasars at high 
redshifts corresponding to the peak epoch of black-hole accretion 
activity, where the host galaxy light in the NIR is sub-dominant. 
The red colours of these quasars in the NIR, can then be attributed 
to dust-reddening rather than galaxy starlight. A WISE colour-cut 
of [Wl — W2)> 0.85 is then applied to separate quasars from 
stars. Candidates are required to be detected at a signal-to-noise- 
ratio (S/N) of >5 in both the WISE W1 and W2 bands. As shown 
in B13, this selection is extremely effective in isolating rare popu¬ 
lations of heavily reddened quasars at high redshifts. 

Our quasar sample encompasses several different survey re¬ 
gions. Below we describe the candidate selection criteria used in 
each of these regions as well as ancillary multi-wavelength data 
that were used to prioritise targets for spectroscopy. After applying 
the colour and morphology selections in each of our survey areas, 
all candidates were visually inspected in order to remove spurious 
sources. The most common such spurious sources are instances of 
close neighbours in the VHS data with separations of <3 arcsec 
and where the WISE photometry appears to be blended, making the 
WISE colours unreliable. The survey regions, selection criteria and 
number of candidates are summarised in Table[T] 

• The VHS-DES StripeSl Region: This region (20h<RA<0h, 
—2° < DEC < 2°) overlaps the SDSS Stripe82 multi-epoch sur¬ 
vey where the coadded data rea ches ~2 mag fainter than the main 
SDSS survey dAnnis et al.l201lh . As discussed in B12, we used the 
deep coadd SDSS Stripe82 photometry to select only the subset of 
red (J — K) candidates that are also faint and red at optical wave¬ 
lengths. Candidates were therefore required to satisfy the additional 
selection criteria: iab > 20.5 and {i — K)ab > 2.5. 

Many of our infrared-selected candidates that failed this cut and 
are blue and bright in the optical, possess SDSS or BOSS spec¬ 
tra. Later, in Section \22\ we show the SDSS spectra to demon¬ 


strate which object populations are eliminated by the optical cuts 
and to determine what the contaminants may be in regions where 
we do not yet have optical coverage. Our selection produced 10 
candidates down to a flux limit of Kab < 18.8 over 129 deg^, af¬ 
ter visual inspection to remove spurious and blended sources. Two 
of the 10, ULASJ2150-0055 and ULASJ2224-0015, were ob¬ 
served as part of the sample described in B12. One more candidate 
(SDSSJ211805.27-1-010344.7) has an optical spectrum from SDSS 
where it is classified as a 2 ~ 5 quasar, although with low confi¬ 
dence. The optical spectrum does not show any obvious emission 
lines. The remaining seven sources were added to our spectroscopic 
target list. 

• The VHS-DES SPT Region: In this region (20h<RA<0h, 
—60° < DEC < —45°) we identify a total of 29 candidates 
down to a flux limit of Kab < 18.8 over 498 deg^. In addi¬ 
tion, over the 83 deg^ area of the ‘SPT Deep Field’ (centred at 
23h<RA<00h and —60° < DEC < —50°), the selection was ex¬ 
tended to TTas < 19.3, giving four additional candidates. The SPT 
Deep Field benefits from multi-wavelength coverage from Spitzer 
dAshbv et alj 1201^ . Herschel and XMM-Newton, which can be 
used to investigate the multi-wavelength properties of our quasars. 

• The VHS-ATIAS SGC Region: In this region covered by VHS 
(20h<RA<0h, —20° < DEC < —2°) we only selected candi¬ 
dates over the ~488 deg^ area with optical photometry from the 
SDSS. The VHS-ATLAS regions have shallowe r J-band coverage 
compared to VHS-DES ( [McMahon et alj[20131) and the selection 
of candidates with extremely red (J — K) colours is prone to con¬ 
tamination from larger numbers of spurious sources. However, un¬ 
like the VHS-DES region, VHS-ATLAS also benefits from Y -band 
coverage, which we used to impose an additional selection crite¬ 
rion: Y > J. Three candidates were selected and verified to be 
faint and red at optical wavelengths based on the SDSS data. 

• The VHS-ATLAS NGC Region: This region 

(13h40<RA<16h20 and -20° < DEC < -2°) encom¬ 

passes the area used for target selection in our pilot survey 
(B13). Nine bright targets with Kab < 18.4 were selected from 
421 deg^ Two of which, VHSJ1350-0503 and VHSJ1409-0830, 
were observed in B13. The remaining seven candidates were added 
to our spectroscopic target list. 

• The UKIDSS-LAS Data Release 10 [H — K) > 1.4 Sam¬ 
ple: Finally, using the UKIDSS-LAS DRIO, covering a total area 
of 3131 deg^ at DEC <20° (observable from the Very Large Tele¬ 
scope), we also constructed a candidate list of the reddest quasars 
which satisfy the additional colour selection criterion of {H — 
K) > 1.4 and have K < 18.9. These would correspond to 
the dustiest hyperluminous quasars at these high redshifts with 
E{B — U) > 1.2 and the objective wa s to find more quasars like 
ULASJ 1234-(-0907 ( [Banerii et al.[[2014l) . Three of these extremely 
red candidates were visible at the time of observations and were 
also added to our spectroscopic target list. 


Our final candidate list comprised a total of 53 spectroscopic 
targets, excluding candidates already observed by B12 and B13. As 
described below, 48 of these targets were successfully observed, 
with 24 classified as broad-line quasars. The candidate list includes 
a complete sample of 41 objects, with Kab < 18.9, over a re¬ 
gion of 1115deg^, selected to provide quantitative constraints on 
the space density of the reddened quasar population. 
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Table 1. Summary of survey regions and selection criteria for reddened quasars presented in this work plus B12 and B13. Several survey regions are spectro¬ 
scopically incomplete. Column 4 gives the number of candidates observed, with the total number of candidates and the percentage completness in brackets. 
Column 5 gives the corresponding information for the number of candidates yielding spectroscopic redshifts, with the total number of spectroscopic targets 
and the percentage redshift yield in brackets. 


Survey Region 

Selection Criterion 

Area/deg^ 

Nobs (Ncand) 

Nw/z (Nobs) 

K < 18.4 Bright Sample 

VHS-DES SPT“ 

{J-K)> 1.6, {Wl - W2) > 0.85 

498 

12(12,100%) 

7(12,58%) 

VHS-ATLAS NGC“’“ 

(J -K)> 1.6, (Wl - W2) > 0.85, Y > J 

421 

5(9,56%) 

3(5,60%) 

VHS-ATLAS SGC“ 

(J -K)> 1.6, (Wl - W2) > 0.85, Y > J 

488 

1(1,100%) 

0(1,0%) 

VHS-DES S82-W“'*' 

(J-K)> 1.6, (Wl - W2) > 0.85 

129 

5(5,100%) 

1(5,20%) 


i > 20.5, (i-K)> 2.5 




UKIDSS-LAS S82-E*' 

(J -K)> 1.6, (Wl - W2) > 0.85 

116 

5(7,71%) 

5(5,100%) 


i > 20.5, (i - K) > 2.5 




UKIDSS-LAS DR10“’‘’ 

(J -K)> 1.6, (Wl - W2) > 0.85, (H -K)> 1.4 

3131 

6(5,83%) 

3(5,60%) 

18.4 < K < 18.9 Intermediate Sample 

VHS-DES SPT“ 

(J -K)> 1.6, (Wl - W2) > 0.85 

498 

17(17,100%) 

9(17,53%) 

VHS-DES S82-W“4 

(J -K)> 1.6, (Wl - W2) > 0.85 

129 

4(4,100%) 

2(4,50%) 


i > 20.5, (i-K)> 2.5 




VHS-ATLAS SGC“ 

(J -K)> 1.6, (Wl - W2) > 0.85, Y > J 

488 

2(2,100%) 

2(2,100%) 

18.9 < K < 19.3 Faint Sample 

VHS-DES SPT Deep“ 

(J -K)> 1.6, (Wl - W2) > 0.85 

83 

4(4,100%) 

1(4,25%) 
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2.2 Optically Bright Reddened Quasar Candidates with 
SDSS Spectroscopy 

Before describing our NIR spectroscopic observations, we examine 
the SDSS/BOSS optical spectra for the optically-bright objects in 
the VHS-DES Stripe82. There are five reddened quasar candidates, 
based on their (J — K) and {W1 — W2) colours, which have 
i < 20.5 and/or (i — K) < 2.5 and were therefore removed from 
our final candidate list. Spectra for these sources are shown in Fig. 
[T] Three objects are at 2 : < 0.6; two appear to be morphologically 
compact star-forming galaxies while one is a low-redshift quasar. 

The two remaining optically bright candidates with SDSS 
spectra are both quasars at higher redshifts. ULASJ2219-I-0036 
(SDSS J221930.42-I-003626.3), a quasar at 2 = 1.196, was also 
identified in B12, where NIR spectroscopy was presented. No emis¬ 
sion lines were detected, confirming that our NIR spectroscopy 
can fail to identify reddened quasars at low redshifts. SDSS 
J220325.00-004002.8, a quasar at 2 = 2.565, possesses a spec¬ 
trum with a blue UV continuum and is therefore not a heavily dust- 
reddened quasar. We note that both these high redshift quasars have 
spectroscopy from the BOSS survey but not from the SDSS Data 
Release 7. _ 

As demonstrated in IRoss et ^ ( l2014h . red quasars selected 
from the SDSS/BOSS surveys constitute a heterogenous popula¬ 
tion including reddened Type 1 quasars. Type 2 quasars. Broad Ab¬ 
sorption Line quasars, as well as sources with extreme emission 
line equivalent widths. Each of these populations is interesting in 
its own right but in this work, we are specifically interested in in¬ 
vestigating the properties of the reddened Type 1 quasars and our 
target selection is designed to isolate this subset of red objects. 


3 SPECTROSCOPIC FOLLOW-UP 

Follow-up of our reddened quasar spectroscopic targets was con¬ 
ducted in Visitor Mode over three full nights in 2013 July using the 
SINFONI spectrograph on the Very Large Telescope (VLT). The 
seeing ranged from 0.5 to l.Oarcsec with a median of 0.78arcsec. 


The quasar candidates were expected to cover a broad range in 
redshift, similar to the samples In B12 & B13, and use of the 
SINFONI H + K filters was therefore requested. However, due 
to technical difficulties, only the SINFONI R=4000 iT-band grat¬ 
ing was available. The 7T-filter has a wavelength range of 1.9- 
2.5/rm, corresponding to a redshift range of 2.1<2<2.7 for Ha 
and redshifts 3.1<2<4.1 for H/3 and [OIII]AA4959,5007. At red¬ 
shifts 0.5 < 2 < 2.1 however, no strong emission lines are visible 
in the iT-band so this represents a redshift desert for our survey. 
ULASJ2219-1-0036 (see Section |2^ is an example of a candidate 
in the redshift desert. 

The dispersion is 2.5A per pixel. Observations were carried 
out in noAO mode, utilising the largest field of view of the instru¬ 
ment, which is 8x8arcsec. The targets were offset by ±1.5 arcsec 
in RA and DEC from the centre of the IFU for sky-subtraction pur¬ 
poses. Apart from the choice of filter, the observational setup is 
identical to that in B12. The total integration time was adjusted de¬ 
pending on the TT-band magnitude of each source and the prevail¬ 
ing conditions at the time of observation. Fainter targets therefore 
have much longer exposure times, chosen such that any broad-lines 
present should be detected at high S/N in the binned spectra. 

All data were reduced using standard ESO pipeline recipes 
using the gasgano data reduction package. Data reduction steps 
included dark subtraction, non-linearity corrections, flat-fielding, 
sky subtraction, extraction and wavelength calibration. Relative 
flux calibration was performed using telluric standards observed 
at similar airmass to the target. Typically, four to five telluric stan¬ 
dards were observed through the night. Due to the absence of spec- 
trophotometric standards in the NIR, absolute flux calibration was 
achieved by normalising the spectra using the JT-band magnitudes 
of the targets. 

From our final list of 53 spectroscopic targets, 48 objects 
were successfully observed with SINFONI over the three nights. 
Twenty-four show broad-lines in the NIR 7T-band spectrum and are 
therefore spectroscopically confirmed to be quasars, with redshifts 
assigned assuming that the emission line is Ha. We note that the 
observed equivalent widths of these lines are inconsistent with cor¬ 
responding to other Balmer lines (e.g. Paa, Pa/3) at lower redshifts. 
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SDSS J220325.00-004002.8 


z=2.565 BOSS QSO 



SDSS J204837.25-002437.2 

z=0.433 BOSS Faint QSO Target 



SDSS J210050.13-005752.5 


z=0.25 BOSS Faint QSO Target 



Rest-Wavelength/A 

SDSS J205601.68-001613.3 


z=0.521 BOSS QSO 



SDSS J221930.42+003626.3 


z=l.196 BOSS QSO 



Figure 1. SDSS spectra of optically bright reddened quasar candidates (black solid) along with the best model-fit used for the redshift determinations (dashed 
red). Note, all of these were quasar targets within BOSS. 


which are expected to be at least an order of magnitude weaker 
dOlikman et al.l200^ . The spectroscopically confirmed quasars to¬ 
gether with their coordinates, iT-band magnitudes, redshifts and 
exposure times are listed in Table One dimensional spectra, ex¬ 
tracted over regions maximising the S/N of the Ha emission line, 
are presented for all objects in Fig. lAll 


4 DEMOGRAPHICS OF THE REDDENED QUASAR 
POPULATION 

4.1 Spectral Energy Distributions, Dust Extinctions & 
Bolometric Luminosities 


© 0000 RAS, MNRAS 000, 000-000 


To provide estimates of dust extinctions and luminosities we fit 
the multi-wavelength photometric data for our confirmed reddened 
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Table 2. Summary of Spectroscopically Confirmed Reddened Quasars. 


Name 

RA 

DEC 

Kab 

Redshift 

Exposure time (s) 

VHS-DES SPT Region 

VHSJ2024-5623 

306.1074 

-56.3898 

18.76 

2.282 

1800 

VHSJ2028-4631 

307.0083 

-46.5325 

18.87 

2.464 

1800 

VHSJ2028-5740 

307.2092 

-57.6681 

17.25 

2.121 

800 

VHSJ2048-4644 

312.0435 

-46.7387 

18.78 

2.182 

1600 

VHSJ2100-5820 

315.1403 

-58.3354 

18.47 

2.360 

1600 

VHSJ2101-5943 

315.3311 

-59.7291 

16.68 

2.313 

320 

VHSJ2115-5913 

318.8818 

-59.2188 

17.63 

2.115 

800 

VHSJ2130-4930 

322.7490 

-49.5032 

18.47 

2.448 

1600 

VHSJ2141-4816 

325.3530 

-48.2830 

18.54 

2.655 

1600 

VHSJ2212-4624 

333.0796 

-46.4101 

18.72 

2.141 

1800 

VHSJ2220-5618 

335.1398 

-56.3107 

16.72 

2.220 

320 

VHSJ2227-5203 

336.9491 

-52.0582 

18.80 

2.656 

2200 

VHSJ2235-5750 

338.9331 

-57.8372 

17.95 

2.246 

1000 

VHSJ2256-4800 

344.1444 

-48.0088 

17.80 

2.250 

1000 

VHSJ2257-4700 

344.2589 

-47.0157 

18.51 

2.156 

1600 

VHSJ2306-5447 

346.5011 

-54.7881 

18.17 

2.372 

1000 

VHSJ2332-5240 

353.0387 

-52.6780 

19.16 

2.450 

2400 

VHS-ATLAS SGC and NGC Regions 

VHSJ1556-0835 

239.1571 

-8.5952 

18.20 

2.188 

1200 

VHSJ2143-0643 

325.8926 

-6.7206 

18.48 

2.383 

1600 

VHSJ2144-0523 

326.2394 

-5.3881 

18.71 

2.152 

1600 

VHS-DES Stripe82 Region 

VHSJ2109-0026 

317.3630 

-0.4497 

18.59 

2.344 

1600 

VHSJ2355-0011 

358.9394 

-0.1893 

18.42 

2.531 

1600 

UKIDSS-LAS DRIO (H - K) > 1.4 

ULASJ0123-I-1525 

20.8022 

±15.4230 

18.59 

2.629 

1200 

ULASJ2315-I-0143 

348.9843 

±1.7307 

18.38 

2.560 

400 


quasa rs using the quasar SED models of iMaddox et al.l 1 I 2 OQ 8 L 
l2012h . In the VHS-DES SPT region the data consists of photometry 
in JHK and WISE W1 and W2. For the Stripe82, VHS-ATLAS 
and UKIDSS regions, optical photometry (from the Stripe82 coadd 
survey and SDSS) is also available and has been utilised in the SED 
fitting. The fitting procedure is described in B12 and accounts for 
the effect of the observed Ha equivalent width on the (J — K) 
colours as follows. The equivalent width of the Ha emission line 
in our base SED model has been derived by fitting to the observed 
colours (ugrizY JHK) of luminous, unobscured, quasars in SDSS 
DR7 that possess SDSS and UKIDSS photometry. The quasar 
SED model (Maddox et al. 2012) takes into account luminosity- 
depen dent changes in the Ha equivalent width due to the Baldwin 
effect (lBaldwiiilll977h . The base quasar-SED model has already 
been shown to provide a very good match to the observed colours 
of luminous SDSS DR7 quasars over the redshift range 0.2<z<4.0 
(Maddox et al. 2012). The default rest-frame equivalent width of 
the Ha line in our z=2.3 model is ~250 A. 

For each of our reddened quasars, we estimate the rest-frame 
Ha equivalent widths directly from the spectra. All equivalent 
widths are calculated by integrating the line within ±3 ct, where a 
is derived from a single Gaussian fit to the line profile. The contin¬ 
uum is defined beyond fhis 3a region out to ±10000 km/s from the 
Ha line centroid. The exact details of the equivalent width calcula¬ 
tions are not important as measurements are are made consistently 
for both the SED model and the observed spectra. The rest-frame 
Ha equivalent widths for each of our reddened quasars are given in 
Table Al and the median value is ~270 A, which is consistent with 
the value for the unobscured quasar population. There is therefore 
no evidence that the reddened quasars have unusual Ha equivalent 
widths relative to their unobscured counterparts. 


There are, however, a few examples of reddened quasars with 
Ha equivalent widths of up to ~500A. For comparison, the K- 
band filter in the VISTA survey has an effective width of 3090 A 
or 940 A in the rest-frame at z=2.3. When deriving E{B — V)- 
estimates, it is important that we account for the effect of the Ha 
line on the continuum iT-band flux. We therefore proceed as fol¬ 
lows. The Ha line in the SED model is scaled to have the same rest- 
frame equivalent width present in each individual reddened quasar. 
As is evident from Table Al the scaling-factor is close to unity for 
the majority of quasars. The E{B — V) values are then derived by 
fitting the scaled SED to all the observed photometry available for 
each quasar. In this way, the specific contribution to the iT-band 
flux by the Ha emission line in each quasar is taken into account 
and the E{B — V) estimates are independent of the emission-line 
strength. 

Our model also includes a luminosity dependent contribution 
from a host galaxy. Even allowing for the effects of extinction on 
the quasar SED, at the redshifts and TT-band magnitudes of our 
sample, the NIR-colours of the majority of the objects are domi¬ 
nated by the quasar SED. Recall that a key element in the sample 
selection is that the candidates must be unresolved in the TT-band 
and, by selection, the host-galaxy contribution cannot be large. The 
form of the host-galaxy SED and exact contribution to the observed 
NIR-colours are thus not significant in the context of reproducing 
the JHKW1W2 photometry. 

The strongest constraints on the reddening for each quasar 
predominantly come from the J — K colour, which probes rest- 
frame wavelengths of ~3400-6400A, below the inflexion in the 
quasar SEDs at ~10000A and above the UV-portion of the spec¬ 
trum where different extinction curves exhibit quite different forms. 
The derived values of E{B — V) are thus completely insensitive to 
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which extinction curve (e.g. curves appropriate for the Milky Way, 
Large or Small Magellanic Clouds) is incorporated in the model 
fits. Given the photometric uncertainties, the estimates of E{B—V) 
are accurate to ~ ±0.1 mag. Extinction values are summarised in 
Tablel^and the median £^(5 — 1/) is found to be 0.8 mag, consistent 
with the values reported in B12 and B13. 

As demonstrated in iBanerii et al.l ( l2014h . these dust extinction 
estimates derived from SED-fitting to the broadband photometry, 
are in excellent agreement with those from independent X-ray ob¬ 
servations, even for the reddest and therefore dustiest quasars. The 
observed quasar SEDs are de-reddene d using th e E{B — V) values, 
assuming an SMC-like extinction law (|PeJ^92h . and the optical lu¬ 
minosity at rest-frame 5100A is calculatecQ. The observed and de- 
reddened i-band quasar magnitudes, together with the de-reddened 
f-band absolute magnitudes and bolometric luminosities are pre¬ 
sented in Table [3 Note, these z-band absolute magnitudes are re¬ 
lated to the SDSS magnitudes via Mi=Miy^^ 2 ] ± 0.596. 

The total luminosities are estimated using a bolometr ic correction 
factor of 7 applied to the rest-frame 5 lOOA luminosity jNetzer et al.l 

I2OO7IR 


where L 5100 is estimated using the best-fit model SEDs of each 
quasar. The Ha rest-frame full-width-half-maximum (FWHM) is 
calculated from either a single- or a double-Gaussian fit to the line 
profile in the velocity range ±10000kms“^. Double Gaussians 
are only used if they produce a statistically significant improve¬ 
ment to the line-profile fit corresponding to a lower mean-fit error. 
All Gaussian components are constrained to possess FWHM in the 
range 1000-10000kms“^, ensuring they trace broad-line region 
gas. Errors on the black-hole masses are derived by propagating the 
mean fit error on the FWHM measurements. Contributions to the 
Ha emission-line profile from any narrow component present are 
small and we verified that including such a narrow Gaussian com¬ 
ponent did not alter significantly the FWHM of the broad single- or 
double-Gaussian components. Parametrization of any narrow Ha 
components, which could potentially be tracing star formation in 
the quasar host galaxy, is discussed in Alaghband-Zadeh et al (in 
preparation). 

The Gaussian line fits are included in Figure lA^ while the 
FWHM values and derived black-hole masses are tabulated in Table 
lAll and Table[3respectively. 


4.2 Black-Hole Masses 


4.3 Comparison to Unobscured Quasars from SDSS 


Single-epoch quasar black-hole masses are normally derived based 
on a virial estimator, employing the widths of broad emission lines 
and assuming a calibration between the size of the broad-line re¬ 
gion and the quasar luminosity. The calibration has bee n derived us¬ 
ing reverber ation mapping techniques at low redshift dKasni et al.l 
I 2 OOOL [ 200 ^ and for relatively low-luminosity AGN; the validity 
of extrapolating the trends to high redshift and higher luminosities, 
still remains uncertain. In addition, virial black-hole mass estimates 
can be sensitive to the choice of broad line used, and the details of 
the line-fitting. A further complication can result if there are sig¬ 
nificant inflows and outflows affecting the gas close to the accret¬ 
ing black-hole that can artificially broaden the observed linewidths. 
If reddened quasars do correspond to the AGN-feedback phase in 
galaxy formation, such inflows and outflows would be more pre¬ 
dominant in our sample and indeed several of the Ha line profiles 
shown in Figure lA2l clearly show evidence for very broad wings 
that could result from such flows. 

Despite these caveats, single-epoch black-hole masses are 
still widely used and there is general consiste ncy between mass 
estimates obtained fr om different broad-lines dShen et al.ll201ll : 
iMatsuoka et al.ll^013h . We therefore provide such mass estimates 
for our redd ened quasar sample. Spe c ificall y, we employ the mass 
estimator of IVestergaard & PetersonI (|2006l), using the ca libration 
between Ha and H/3 linewidths from iGreene & Ha ( l2005h : 


FWHMh/ 3 = (1.07±0.07)xl0® / FWHMna 

V lO^kms 


kms 


( 1 ) 


MBH/Me = 10-- ^ (Jr " (2) 

VlOOOkms V 


^ As with the determination of the E{B — F)-values, the derived opti¬ 
cal luminosities are insensitive to which conventional extinction curve is 
employed. 

^ A holo metric correction factor of closer to 9 has also been used in the 
literature IShen et alj|201lh and would increase our estimated bolometric 
luminosities by ~0.1 dex. 


The black-hole masses and bolometric luminosities can be used 
to derive Eddington ratios (L/L-Edd) which are typically >0.1, 
as was found in B12 and B13. The reddened quasars therefore 
appear to be massive black holes, accreting at a high rate, re¬ 
sulting in their very high luminosities. Large numbers of unob¬ 
scured, luminous and highly accreting black-holes are known, 
e.g.fr om SDSS jSchneider et alj 1 20071) and BOSS jParis et ^ 
l2014h . and the luminosity function and space density of the 
obscured and unobscured populations are of potential interest. 
We define a sample of 2081 ‘unobscured’ quasars from SDSS 
DR7, with redshifts 2.1< z <2.7 and magnitudes K <18.9 
as w ell as black - hole m ass and bolometric luminosity estimates 
from IShen et ^ l l201lh . The combined sample of spectroscopi¬ 
cally confirmed reddened quasars from B12, B13 and this paper, 
over the same redshift range and brighter than K <18.9 repre¬ 
sents the ‘reddened’-quasar population. The distribution of black- 
hole masses and luminosities for the two samples are shown in 
Fig.lH The mean black-hole masses and bolometric luminosities 
are logio(MBH/M0)=9.7±O.4, logic(Lboi/erg s“^)=47.1±0.4 and 
logio(MBH/M0)=9.3±LL logic(Lboi/erg s~^)=47.0±0.2 for the 
reddened and unobscured samples respectively. 

The distributions of black-hole mass and luminosity overlap 
significantly (Fig. [2j, but there is an indication that the reddened 
quasars possess a tail out to higher bolometric luminosities and 
black-hole masses. Several selection effects and systematic errors 
could be responsible for this observed tail. The black hole masses, 
derived from UV-emission lines for objects in SDSS and Ha for 
the reddened population, may not be on exactly the same scale. 
Progress is being made via direct comparison of emission-line 
prope rties in quasar samples jMatsuoka et al.|[2013l : lTilton & Shulll 
l2013h but the uncertainty in black hole mass estimates from in - 
dividual emission lines remains significant (e.g. lShen et al.ll2oTlh . 
More importantly, inflows and outflows are likely to be more com¬ 
mon in the reddened population, thought to be caught in the AGN 
radiative feedback phase. The presence of such flows could add 
significant additional spread to the observed linewidths, increasing 
black-hole masses by factors of a few (B12). The black-hole masses 
should therefore be interpreted with caution. 

The bolometric luminosities of the samples may also differ 
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Figure 2. Distribution of bolometric luminosity (left) and black-hole mass (right) for the red and unobscured quasar populations (see descriptions in Sec¬ 
tion 4.3). There is significant overlap between the samples in both luminosity and mass, although there is a suggestion that the red population includes a tail of 
quasars with very high bolometric luminosities and black-hole masses. 


systematically. Luminosities for the unobscured quasar sample are 
calculated using the 1350A rest-frame luminosity and assume a dif¬ 
ferent model quasar SED to ours. Most likely, however, is a selec¬ 
tion against the most luminous quasars in the SDSS resulting from 
the bright z-band flux-limit of 15.0 mag (c.f. the values for the red¬ 
dened quasars in Col. 4 of Table [3]l. 

Overall, notwithstanding the caveats summarised above, we 
can conclude that the reddened quasar sample consists of a popula¬ 
tion of very high mass black-holes, comparable to the most massive 
black-holes included in the SDSS quasar sample. 


5 INFRARED SEDS EROM WISE 

Comparing the multi-wavelength properties of the reddened quasar 
sample to other well-studied galaxy and AGN populations at simi¬ 
lar redshifts helps in understanding the origin of the sample. The 
availability of the WISE photometry, probing rest-frame wave¬ 
lengths out to >3 fim for objects detected in the VF3-band, provides 
information for a key portion of quasar SEDs. An unobscured com¬ 
parison sample is defined by selecting 1604 SDSS DR7 and DRIO 
quasars with 2.1 < 2 < 2.7, Kab < 18.9 and detections in the 
WISE VL123-bands - hereafter the ‘SDSS sample'. 

The reddened quasar sample was selected using a (W1 — 
W2) > 0.85 colour cut as detailed in Sectionl^and are therefore, 
by definition, red in the bluer WISE bandfl We show the WISE 
colours of our spectroscopically c onfirmed quasars in Eig.[^on the 
colour-colour locus introduced bv lWright et al.l ( l2010h . All the ob¬ 
jects are detected at > 3a in the W3 12 /rm band and the majority 
have W lVE2TE3-colours placing them in the QSO/Seyfert locus of 


® All WISE photometry presented in this section are taken from the new, 
improved WISE data reductions available through the AUWISE data release. 
However, the WISE main survey catalogue was used for the target selec¬ 
tion as the AUWISE reductions were not available at the time. We find that 
the WISE and AUWISE fluxes are generally consistent apart from a single 
quasar: VHSJ2024-5623. The quasar satisfied the (W1 — W2) > 0.85 
colour cut in the original WISE data release with W1 = 16.52 and 
W2 = 15.49 but has revised magnitudes in AUWISE of W1 = 16.17 
and W2 = 15.44. 


IWright et alj ( l2010l) . Just a few of the reddest sources extend into 
the Obscured AGN locus. 

Eig.|3]suggests that, crudely, the reddened quasars are indistin¬ 
guishable from unobscured quasars in terms of their WISE colours. 
We can test this quantitatively by fitting a power-law to the ob¬ 
served WISE fluxes in fhe Wl, W2 and W3 bands of the form 
XFx oc . Emission from hot dust is expected to dominate the 

SEDs over the 1^ /zm rest-frame wavelength range probed by the 
WISE passbands. We find a mean value of /3nir=0.69 in our sam- 
plt0 with a standard deviation of 0.29. The mean NIR power-law 
inde x for the SDSS sam ple i s 0.50±0.26, cons istent with the results 
from f^ng et aD ( l2013h and iHao et alj ( 1201 ll) . 

The SDSS and reddened-quasar power-law indices agree 
within Ict of the population deviation and the steeper mean value 
for the reddened quasars is consistent with the most obscured of 
our quasars showing some extinction at rest-frame wavelengths of 
l-2^m, which steepens the slope of the SED. Eig.|4]shows the dis¬ 
tribution of 13 for both populations. The high /Jnir tail is due to the 
reddest quasars in our sample, with E{B — V^) > 0.8, as expected. 

Several recent studies have focused on the subset of high 
redshift AGN tha t are extremely bright in the W4: 22 /zm band. 
IRoss et aT] (1201 4l) studied a sample of SDSS and BOSS spectro¬ 
scopically confirmed quasars wifh W4 <= 8 mag and the WISE 
consortium have selected hyperluminous infrared galaxies, with 
W4 < 7.7 mag, that are extremely red in their {W1 — W2) colours 
dEisenhardt et alj2012l : IWu et al.l2oT^ . Our selection for reddened 
quasars does not use the 22 /zm fluxes, allhough, by selecting ob¬ 
jects that are bright in the JT-band, we are by definition not sensi¬ 
tive to the WISE ‘W12’ drop-outs. Twenty-five of our 38 reddened 
quasars are detected at > 3cr in the W4 band. Of these, 14 have 
W4 <= 8 mag and 11 have W4 <= 7.7mag, 37 and 29 per cent 
of the sample respectively. The 22 /zm-detected objects are gener¬ 
ally the most bolometrically luminous quasars in our sample, as 
expected. 

Eor objects detected in the WISE VF4-band, rest-frame lu¬ 
minosities at 6/zm can be estimated using the W3 and W4 


The convention for /3 nir employed here is in wavelength rather than 
frequency space compared to our previous work (B13). The two slopes are 
related by /3nir=/3nir(B13)-h1. 
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Table 3. Best-fit dust extinction, reddened and de-reddened optical i-band magnitudes, absolute i-band magnitudes, bolometric luminosities, black-hole masses 
and Eddington ratios of the reddened quasars. 


Name 

E{B-V) 

:red 

^AR 

^clerecl 

^AR 

Mi 

logio(Lboi/erg s 

logic (Mbh/M©) 

lEdd 

dm/dt (M 0 yr 

VHS-DES SPT Region 

VHSJ2024-5623 

0.6 

22.9 

17.8 

-27.37 

46.7 

9.8±0.1 

0.06 

9 

VHSJ2028-463I 

0.6 

22.9 

17.9 

-27.49 

46.7 

9.2±0.2 

0.2 

9 

VHSJ2028-5740 

1.2 

24.3 

14.6 

-30.18 

47.9 

lO.liO.l 

0.4 

127 

VHSJ2048-4644 

0.7 

23.2 

17.4 

-27.65 

46.8 

9.2±0.1 

0.3 

11 

VHSJ2100-5820 

0.8 

23.6 

16.5 

-28.49 

47.1 

9.1 ±0.1 

0.8 

24 

VHSJ2101-5943 

0.8 

21.8 

14.9 

-30.03 

47.8 

10.5±0.1 

0.1 

106 

VHSJ2115-5913 

1.0 

23.5 

15.4 

-29.48 

47.5 

9.3±0.1 

1.4 

61 

VHSJ2130-4930 

0.9 

24.7 

16.5 

-28.75 

47.2 

9.5±0.2 

0.4 

31 

VHSJ2141-4816 

0.8 

23.8 

16.1 

-29.10 

47.3 

9.6±0.2 

0.4 

40 

VHSJ2212-4624 

0.8 

23.7 

17.2 

-27.64 

46.9 

9.9±0.1 

0.07 

14 

VHSJ2220-5618 

0.8 

21.6 

14.9 

-30.09 

47.8 

9.9±0.1 

0.5 

106 

VHSJ2227-5203 

0.9 

25.0 

16.3 

-29.15 

47.4 

10 . 0 ± 0.1 

0.2 

41 

VHSJ2235-5750 

0.6 

21.9 

16.9 

-28.28 

47.1 

10 . 1 ± 0.1 

0.07 

20 

VHSJ2256-4800 

0.6 

21.7 

16.7 

-28.49 

47.1 

10 . 1 ± 0.1 

0.08 

24 

VHSJ2257-4700 

0.7 

23.0 

17.3 

-27.67 

46.9 

9.5±0.1 

0.2 

13 

VHSJ2306-5447 

0.7 

23.2 

17.1 

-28.15 

47.0 

10 . 0 ± 0.1 

0.08 

18 

VHSJ2332-5240 

0.6 

23.4 

18.0 

-27.39 

46.7 

9.5±0.1 

0.1 

8 

VHS-ATL SGC and NGC Regions 

VHSJ1556-0835 

0.7 

22.7 

16.9 

-28.15 

47.0 

8 . 6 ± 0.1 

1.9 

19 

VHSJ2143-0643 

0.8 

23.8 

16.7 

-28.33 

47.1 

10 . 0 ± 0.1 

0.09 

22 

VHSJ2144-0523 

0.6 

22.6 

17.7 

-27.33 

46.7 

9.9±0.1 

0.04 

9 

VHS-DES Stripe82 Region 

VHSJ2109-0026 

0.7 

23.3 

17.2 

-27.95 

46.9 

9.8±0.1 

0.1 

15 

VHSJ2355-0011 

0.7 

23.2 

16.6 

-28.67 

47.2 

10 . 1 ± 0.1 

0.09 

27 

UKIDSS-LAS DRIO (H - K) > 1.4 

ULASJO123-1-1525 

1.3 

27.7 

15.4 

-30.12 

47.8 

9.7±0.2 

0.9 

114 

ULASJ2315-1-0143 

1.1 

26.2 

16.1 

-29.43 

47.5 

10 . 1 ± 0.2 

0.2 

57 


Table 4. WISE W3 and W4 magnitudes for all reddened quasars detected at > Scr at 22/im along with the rest-frame 6 /im luminosity. 


Name 

W3 (mags) 

iy4(mags) logio{L 6 iim/erg s 

ULASJO 123-h 1525 

8.69±0.03 

6.33±0.06 

47.8 

ULASJ0144-0114 

10.07±0.05 

8 . 02 ± 0.21 

47.1 

ULASJ0221-0019 

10.46±0.06 

8.11±0.17 

46.9 

ULASJ1234-H0907 

9.17±0.04 

7.28±0.17 

47.4 

VHSJ1350-0503 

10.62±0.07 

8.59±0.26 

46.7 

VHSJ1409-0830 

10.33±0.06 

8.31±0.23 

46.9 

ULASJ1455-H1230 

10.43±0.05 

8.11±0.15 

46.3 

ULASJ1539-H0557 

9.10±0.03 

6.56±0.06 

47.8 

VHSJ2028-5740 

9.15±0.03 

7.01±0.08 

47.3 

VHSJ2048-4644 

10.42±0.08 

7.92±0.18 

46.9 

VHSJ2100-5820 

10.70±0.10 

8.14±0.24 

46.9 

VHSJ2101-5943 

9.25±0.03 

6.84±0.07 

47.4 

VHSJ2115-5913 

9.48±0.04 

7.06±0.08 

47.2 

VHSJ2130-4930 

10.69±0.09 

8.52±0.33 

46.9 

VHSJ2141-4816 

10.46±0.07 

7.68±0.16 

47.3 

ULASJ2200-H0056 

10.23±0.06 

7.80±0.17 

47.2 

VHSJ2212-4624 

11 . 02 ± 0.12 

8.37±0.27 

46.7 

VHSJ2220-5618 

8.77±0.03 

6.37±0.05 

47.6 

VHSJ2227-5203 

10.11 ±0.06 

8 . 02 ± 0.22 

47.2 

VHSJ2235-5750 

9.97±0.05 

7.84±0.17 

47.0 

VHSJ2256-4800 

9.67±0.04 

7.52±0.13 

47.1 

VHSJ2257-4700 

10.77±0.09 

8.43±0.28 

46.7 

ULASJ2315-H0143 

9.57±0.04 

7.03±0.10 

47.5 

VHSJ2332-5240 

11.51 ± 0.20 

8.45±0.31 

46.9 

VHSJ2355-0011 

9.83±0.05 

7.40±0.14 

47.3 
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Figu re 3. Location of red dened quasars in WISE colour-colour plane taken 
from lWright et alj llOlOl) showing that they overlap the QSO/Seyfert locus. 



Figure 4. Distribution of NIR power-law indices, /3 nir, fit in the rest-frame 
wavelength range 1-4/rm both for our reddened quasar sample and SDSS 
unobscured quasars over the same redshift range and down to the same K- 
band flux limit. 


photometry. At the median redshift of z =2.3 for our sample, 
the W3 hand corresponds to 3.6/rm and the W4 hand corre¬ 
sponds to 6.7 /im. The 6 /rm luminosities are calculated via lin¬ 
ear interpolation, using the rest-frame wavelengths correspond¬ 
ing to the W3 and W4 magnitudes, and the results are given in 
Table |4] and illustrated in Fig. |5] The average rest-frame 6 
luminosity is logiofLe^m/erg s“^)=47.1±0.4. For a comparable 
SDSS DR7-I-DR10 sample of 1336 2.1< z < 2.7 quasars 
with K< 18.9 and >3a detections in the FF4 band, the av¬ 
erage luminosity is logiofLe^m/erg s“^)=46.8±0.2. There are 
six SDSS quasars with logiofLe^m/erg s~^)> 47.5 (four from 
SDSS DR7 and two from SDSS DRIO). The brightest of these 
six is the well-kno wn gravitationally-lensed Cloverleaf quasar 
jMagain etal.ll988h . SDSSJ000610.67-I-121501.2 (z = 2.309) and 
SDSSJ155102.79-I-084401.1 (z = 2.520), the two DR IP objects 
were p resented in the extremely red sample studies bv iRoss et al.l 
(l2014f) . By comparison, our reddened quasar sample, selected us¬ 
ing imaging data covering less t han a third of t he total SDSS 
imaging area 14,555 deg^ from IPms et alj|2014h . includes four 
quasars with logiofLe^m/erg s“^)> 47.5, highlighting the in¬ 
creased prevelance of such luminous quasars among the obscured 
sample. Compared to VV75£-selected HyLI RGs where the most lu - 
minous object at 6 /rm, WISEJ1814-b3412j [Eisenhardl^al.l2012h . 
has logiofLe^m/erg s“^)=47.30±0.05 jStem et alj|2014h . the red¬ 
dened quasar sample is also extreme, with seven objects more lu¬ 
minous than WISEJl814-1-3412. 

Einally, we compare the 6 fj,m luminosities to the bolometric 
luminosities in Table and find a typical bolometric correction 
factor of ~2, consistent with a model in which a large fraction of 
the total luminosity is re-radiated at mid infrared wavelengths. This 
bolometric correction from the mid infrared is also broadly consis¬ 
tent with that derived for the most luminous essentia lly unobscured 
quasars in SDSS and WISE jWeedman et al.ll2012h . As the 6 fim 
rest-frame luminosities are inferred directly from the observed pho¬ 
tometry without any assumptions regarding the dust extinction law 
and the form of the quasar SED, this analysis corroborates our 
conclusion in Section IdTl that our reddened quasars correspond 
to some of the most intrinsically luminous quasars known. 



I • SDSSDR7+DR10 
• Red Quasars 
I ■ WISE HyLIRGs 


48.2 

48 

47.8 


Redshift 


Figure 5. Redshift versus rest-frame 6/rm luminosity for our reddened 
quasars compa red to unobscured quasars from SDSS as well as 3 WISE 
HyLIRGs from IStem et alj j2014h . The brightest unobscured quasar on this 
plot is the lensed Cloverleaf quasar. 


6 OBSCURED BLACK-HOLE GROWTH AT THE 
HIGHEST QUASAR LUMINOSITIES 

With a sample of 38 reddened quasars {E{B — V) >= 0.5), of 
which 36 lie at z > 2, now in place, a direct comparison of their 
space density to that of unobscured quasars can be made. As deep 
X-ray and mid-infrared surveys have focussed on studying much 
fainter, dusty, AGN, the sample allows us to address the key ques¬ 
tion of how many very luminous quasars are dust-obscured at red- 
shifts 2 ~ 2 — 3. We begin by considering the selection function 
and completeness of our reddened quasar sample. 


6.1 Sample Completeness 

6.1.1 Morphology Cut 

The sample selection (Section |2j includes a requirement that can¬ 
didates are unresolved in the Tf-band. While eliminating contami¬ 
nation by low-redshift quasars, where host galaxies can contribute 
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significantly to the observed colours (B12), it is necessary to quan¬ 
tify the fraction of high-redshift quasars that may also be ex¬ 
cluded. To assess such incompleteness, we select a sample o f un¬ 
obscured quas ars from the SPS S DR7 l lSchneider et aLlEoiOh and 
SDSS DRlO l lParis et an i2014h quasar catalogues, with redshifts 
2.1 < 2 : < 2.7 and magnitudes K < 19.3. There are a total of 
7871 quasars that satisfy these criteria. The cumulative fraction of 
unobscured quasars classified as unresolved as a function of K- 
band magnitude is shown in Fig. [7] At the bright-end (K < 18), 
almost 90 per cent of the quasars are unresolved, with the fraction 
decreasing to just under 70 per cent at the faint limit. Only four of 
the 48 reddened quasar candidates are in the faintest magnitude bin 
and, using the distribution of candidate iT-magnitudes, we calcu¬ 
late the incompleteness due to the image morphology constraint to 
be 20-25 per cent. 


Redshift 
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6.1.2 (J — K) Colour Cut 

The NIR colour-cut of (J — K) > 1.6 is designed to isolate 
the most highly reddened Type 1 broad-line quasars. Fig. |6] shows 
model quasar (J — K) colours, calculated using our quasar SED 
model (Section |4. lb . over the redshift range of our sample, as a 
function of the extinction parameter E(B — V). Some dependence 
of (J — K) colour on redshift is present; higher redshift quasars 
have bluer {J — K) colours due to the Ha emission line moving 
out of the iT-band. However, there is a tight correspondence be¬ 
tween extinction and colour, with the (J — K) > 1.6 threshold 
predicted to define a quasar sample with E(B — H) > 0.5. 

A potential concern is that quasars with very large Ha equiv¬ 
alent widths and somewhat lower E{B — V) values form a signif¬ 
icant fraction of the sample. The direct measurements of the Ha 
equivalent widths for the reddened quasars (Table IaH , however, 
demonstrate that the presence of such objects is not a factor. The 
highest equivalent width quasar in our sample has an equivalent 
width of only a factor of ~2 larger than the average equivalent 
width in unobscured quasars of similar luminosity. The average Ha 
equivalent widths however are very similar to those seen in unob¬ 
scured quasars, confirming that the reddened quasars we have se¬ 
lected do not have unusual equivalent widths relative to the unred¬ 
dened population. 

The systematic dependence of the threshold E{B — H) as a 
function of redshift to satisfy the (J — K) > 1.6 selection, evi¬ 
dent in Fig. 1^ is directly taken into account when estimating the 
space-density (Section l6.3t . The S/N of the iT-band photometry 
used in the sample selection is high and any apparent increase in 
the ob served number of objects due to Eddington bias dEddinetonl 
Il913h is small. The S/N of the J-band photometry is somewhat 
lower and the E(B — V) selection at a given redshift, resulting 
from the (J — K) > 1.6 colour cut, is thus blurred, with some 
objects lost from the sample and others included due to photomet¬ 
ric errors. This could be a concern if the E{B — V) distribution 
was changing significantly at the boundary of our selection box 
{E{B — H) ~ 0.5), in which case unequal numbers of objects 
may be scattered out of the sample relative to the numbers of ob¬ 
jects that are scattered in. Howeve r, estimates of the fo rm of the 
E{B — V) distribution from SDSS I Hopkins et al.l200^ and near- 
IR selected samples dMaddox et al.ll2012ll strongly sueeest that the 
E{B — V) distribution is not changing significantly and our sam¬ 
ple should thus represent an accurate estimate of number of objects 
with E{B — V) exceeding the limits calculated. 


Figure 6. Observed, reddened (J — K) colour as a function of extinction 
E{B — V) for a set of simulated quasars assuming the base SED model 
described in Section lTTI and an SMC-like extinction law. Points are colour- 
coded by the redshift of the quasar. The dashed horizontal line shows our 
colour selection for the reddened quasars. 


6.1.3 WISE Colour Cut 

We employ the same SDSS sample of unobscured quasars used 
in Section [6. 1.1 1 to determine the incompleteness of the reddened 
quasar sample due to the application of the restriction that WISE 
(Wl — W2) >= 0.85. The WISE colours of the reddened quasar 
sample are themselves affected by the presence of dust. Specifi¬ 
cally, at the median redshift of the sample (z = 2.3), the W 1 and 
W2 magnitudes of a quasar possessing an E{B — V) of 0.5 mag 
are 0.75 and 0.53 mag fainter respectively. The requirement that 
{W1 — W2) >= 0.85 forreddened quasars with 77(73 —H) > 0.5 
therefore corresponds to {W1 — W2) >= 0.63 for the unobscured 
SDSS sample. 

Eigure|7]also shows the cumulative fraction of the unobscured 
SDSS sample with (Wl — W2) >= 0.63, W1s/n > 5 and 
W2g/iv > 5 (dot-dashed line). This fraction is >90 per cent for 
most of the magnitude range but falls to ~85 per cent in the faintest 
magnitude bin, used in the 83 deg^ VHS-DES SPT Deep Field sur¬ 
vey. 

The cumulative fraction of quasars, with redshifts 2.1 < 2 < 
2.7, that satisfy both the (J — K) and (Wl — W2) colour cuts, 
plus the unresolved morphology constraint, i.e. the combined sam¬ 
ple photometric completeness, is also shown in Fig. |7] More than 
70 per cent of quasars satisfy all our selection criteria at K = 18.9 
with the completeness falling to just under 65 per cent atK = 19.3. 


6.2 Redshift Completeness 

In addition to the candidate photometric completeness, the effec¬ 
tiveness of the spectroscopic observations for determining red¬ 
shifts and hence reliable object classifications must be considered. 
Fig. lAll shows that the high equivalent-width, broad, Ha emission 
lines are clearly detected at high S/N, even at the faint 7T-band limit 
(e.g. VHSJ2332-5240) or when the line is at the edge of the wave¬ 
length interval (e.g. VHSJ2227-5203). All Ha lines have S/N> 4. 
The observational strategy (Section[3](, where integration times (Ta¬ 
ble |2j were adjusted to achieve a minimum S/N for each target 
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Figure 7. Cumulative fraction of SDSS DR7 and DRIO quasars at 2.1 < 
2 < 2.7 and K < 19.3, that are also classified as point sources in the K- 
band, as a function of JC-band magnitude (dashed line). The dot-dashed line 
shows the cumulative fraction of SDSS quasars that have (W1 — W2) > 
0.63 (corresponding to (W1 — W2) > 0.85 for E{B — V) > 0.5). The 
solid line shows the cumulative fraction of SDSS high redshift quasars that 
satisfy both the moiphology and WISE colour selection. The solid vertical 
line marks the magnitude limit used for the majority of our reddened quasar 
sample while the dashed vertical line shows the magnitude limit used in the 
83 sq-deg VHS-DES SPT deep field. We can see that >70 per cent of true 
high redshift quasars would satisfy all our selection criteria at K=18.9 and 
this completeness falls to just under 65 per cent at K = 19.3. 


sufficient to allow detection of broad emission lines, thereby en¬ 
sures that the redshift completeness is close to 100 per cent within 
the redshift interval 2.1 < 2 < 2.7. Fig.[^shows the rest-frame 
equivalent width of the Ha line as a function of the quasar K-hand 
magnitude. No significant correlation between the two quantities is 
present and, particularly for the faintest quasars, a range of equiva¬ 
lent widths is evident. There is no evidence to suggest that, for the 
faintest sources, our survey is preferentially detecting quasars with 
large Ha equivalent widths. 

In the space-density analysis that follows, we first consider 
a ‘Bright Sample’, from the VHS-DES SPT, VHS-DES S82-W 
and VHS-ATLAS SGC regions, as summarised in Table Q] These 
three regions cover a total area of 1115 deg^ and all 29 candidates 
with Kab < 18.9 possess spectra. Twenty one of the 29 candi¬ 
dates have redshifts. In addition. We also consider a ‘Faint Sam¬ 
ple’, based on the Kab < 19.3 sample, covering the 83 deg^ SPT 
Deep Field, where, again, all candidates possess spectra. 


6.3 The Reddened Quasar Space Density 

We calculate the V/Knax weighted dSchmidd 1 19681 : 
lAvni & Bahcail Il980h estimate of the reddened quasar space 
density from the 21 reddened quasars in the Bright Sample. Using 
the hest-fit SED for each reddened quasar, we can calculate its 
absolute magnitude as well as the redshift range and therefore 
total volume V, within which it would have satisfied our survey 
selection criteria. Our space density calculations assume the quasar 
SED model from Section 14.11 in order to derive k-corrections 
and therefore absolute magnitudes and bolometric luminosities. 
As our SED model is slightly different from the one used in 


Figure 8. Rest-frame equivalent width of the Ho line in our reddened 
quasars as a function of the if-band magnitude. The horizontal dashed 
line represents the average Ho equivalent widths of unobscured quasars 
at the same redshift and luminosity. The Ho lines have a range of equiva¬ 
lent widths between ~ 0.5 — 2 times the equivalent widths in unobscured 
quasars. Particularly at the faint-end of the sample, a range of Ho equivalent 
widths is evident with no strong coirelation between the A-band magnitude 
and the equivalent width. 


SDSS teoss etalj|20l3l) and in order to allow comparison with 
investigations by other groups, we provide a table with the i-band 
and A-band k-correction from our model in Table m However, 
we note tha t differences in k-c orrections derived from our template 
and that of IRoss et alJ l l2013t) are expected to be small given that 
both tem plates are able t o repr o duce the observed colours of SDSS 
quasars ( iMaddox et alj 1201 2l : IRoss et ah I I 2 OI 3 I) . The maximum 
available survey volume, Vjnax is the co-moving volume within 
2.1 < z < 2.7. The cumulative volume and completeness- 
corrected number density in each absolute magnitude bin is then 
given by: 


N{M < Mi)carr/MpC® 


V max 


Y^Nbin 


1 

c(K) X Vj/V^ax 


(3) 


The error on these cumulative number counts is: 




fsky ^Nbin f 

Vma^ " = ^ (C(K) X Vj/Vmax)2 


(4) 


where fsky denotes the fraction of area on the sky encompassed 
by our survey and c{K) is the A-band magnitude dependent pho¬ 
tometric completeness of each quasar calculated using the com¬ 
pleteness curve in Eig.|3 The resulting cumulative space density is 
shown in Fig. [^together with the equivalent space density for un¬ 
obscured c|uasare_^rived from the SDSS-l-BOSS luminosity func¬ 
tion (LF: IRoss et alj[20l^ over exactly the same redshift and ab¬ 
solute magnitude range as the reddened sample. As our luminosity 
function for reddened quasars extends to brighter absolute mag¬ 
nitudes than the SDSS-l-BOSS quasars, we assume the parametric 
double power-law form of the unobscured quasar luminosity func¬ 
tion and integrate this parametric luminosity function over the ap¬ 
propriate redshift and luminosity ranges to compare to the reddened 
quasars in Fig. We consider two best-fit forms for this para- 
metric luminosity function for unobscured quasars from iRoss et all 
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( l2013h . Firstly, a luminosity and density evolution (LEDE) model 
and secondly a pure luminosity evolution model (PEE). We note 
that the SDSS+BOSS luminosity function is constrained in the red- 
shift range 2.2 < z < 3.5 and only a very modest extrapolation of 
the luminosity function and redshift-evo l ution is necessary to com¬ 
pare with our new data. In IRoss et all (l2013h . a pure luminosity 
evolution (PEE) model has instead been fit to the data al z < 2.2 
but the effect on the cumulative LEDE model based number counts 
of the unobscured quasars plotted in Fig. is small. We also note 
that the systematic differences in k-corrections between our models 
and those used by BOSS, are expected to be smaller than the bin 
size in Fig.|^ 

We find that the LF of the reddened quasars is significantly 
flatter than that of the SDSS-l-BOSS quasars over the entire magni¬ 
tude range. At bright absolute magnitudes (Mi < —29), the red¬ 
dened quasars possess larger space densities than predicted by the 
SDSS-hBOSS LEDE LF, but the more steeply rising SDSS-l-BOSS 
quasar LF means the unobscured quasars become significantly 
more common at magnitudes Mi > —27.5. The much flatter in¬ 
crease in the number of quasars as a function of decreasing lu¬ 
minosity leads to very different interpretations regarding the frac¬ 
tion of the total energy liberated by the Type 1 AGN popula¬ 
tion in an ‘obscured’ phase, compared to previous i nvestigations 
which probed only the mo st luminous quasars (e.g. IWhite et aP 
I 2 OO 3 I : [oiikman et al.ll2007h . Our survey now covers enough dy¬ 
namic range in luminosity to demonstrate that while the numbers 
of obscured (0.5 < E{B — V) < 1.5) quasars are significant at 
the highest luminosities, their space density shows a much shal¬ 
lower rise relative to unobscured quasars at these redshifts. There 
is no evidence from our data that the reddened quasars have signif¬ 
icantly different Ha equivalent width distributions relative to un¬ 
obscured quasars so we are unlikely to be missing a population of 
very low equivalent width objects but with similar levels of redden¬ 
ing as in our sample. As such, these reddened broad-line quasars 
are likely to be sub-dominant relative to the unobscured quasars 
in terms of their contribution to the total energy liberated during 
the luminous quasar phase. For reference, our K =18.9 magnitude 
limit corresponds to Mi ~ —27.5 at z = 2.3 for a quasar with 
E{B — V) = 0.75. This is ~1.3mag brighter than M* at these 
redshifts, where the bulk of the luminosity is produced. A larger 
survey to fainter iT-band limits is required to determine the ex¬ 
act fraction of reddened quasars down to the knee of the intrinsic 
quasar luminosity function and hence quantify what fraction of the 
accretion luminosity for the Type 1 AGN-population is generated 
in an obscured phase. 

Our result, that reddened quasars are more common at bright 
absolute magnitudes, contrasts with observations of intrinsically 
less luminous AGN-populations, where the obscured fraction 
shows the opposit e trend and decrease s with increasing lu minosity 
jUeda et al.ll2003tl . On the other hand. lAssef et alj ( |2014|) have re¬ 
cently studied the space density of the highly obscured population 
of WT5£-selected HyLIRGs, which overlap our sample in terms of 
both luminosity and redshift, but have higher extinctions. They also 
find that the fraction of obscured luminous AGN is significant at 
very high luminosities. 

Our sample is not large enough to merit fitting a parametric 
form for the evolution of the reddened quasar LF as a function 
of redshift. Nevertheless, the differences in the cumulative number 
counts of unobscured quasars seen from the LEDE and PLE LFs 
in this redshift range, can help shed light on the kind of evolution¬ 
ary model that may be appropriate. Starting with the SDSS-l-BOSS 
luminosity functions in Fig. [9] and assuming that (i) the extinction 


is independent of luminosity and (ii) that every quasar possesses 
a reddening between 0.5<E{B — V)<2.0 with reddening values 
drawn randomly according to a flat probability distribution in this 
range, we can estimate how many reddened quasars we would have 
observed in our survey area employing the completeness function 
shown in Fig. |7] These predicted counts, when compared to our 
observed counts, can then help determine whether the underlying 
luminosity distribution of reddened quasars, is indeed the same as 
that of unobscured quasars. Based on the normalisation provided by 
the BOSS LEDE LF, we predict 15±4 reddened quasars would be 
detected in the Bright Sample, which is lower than the 21 quasars 
we have observed. Adopting in stead the E(B — H) d istribution in¬ 
ferred from SDSS quasars by iHopkins et al.l l l2004h rather than a 
flat E{B — V) distribution in the range 0.5<E{B — V)<2.0, re¬ 
sults in 0.35±0.59 quasars predicted in our survey. The dramatic 
difference in the number of objects provides a direct illustration of 
the need to incorporate near-IR photometric selection in order to 
identify quasars with reddenings of E{B — I/)~0.5-L0. Finally, 
adopting the BOSS PLE LF model instead, which has a brighter 
break-luminosity (M*) than the LEDE model, results in a larger 
space density for the most luminous quasars as evident in Fig.|3 As 
a consequence, using the normalisation from the PLE LF-evolution 
model would predict a larger number of reddened quasars in our 
Bright Sample - 34±4. Neither model, together with the assump¬ 
tion that the extinction is independent of luminosity, provides a 
good match to our observed number counts of reddened quasars 
. The SDSS-l-BOSS LEDE LF produces too few reddened quasars 
while the PLE LF produces too many. The results demonstrate that 
reddened quasars with —27.0 < Mi < —30.0 do not possess the 
same intrinsic [unreddened] LF as unobscured quasars unless the 
intrinsic luminosity and the E{B — V) distribution are related. 
Qualitatively, the probability that a quasar possesses a significant 
E{B — V) must be greater for intrinsically more luminous objects 
in order for the SDSS-l-BOSS luminosity distributions to be consis¬ 
tent with our data. The quasar luminosity and reddening are thus 
not independent. A model in which the probability of a quasar hav¬ 
ing high E{B — V) increases with increasing luminosity would be 
viable. However, in the conventional AGN unification model, in¬ 
corporating a receding torus, the dependence is predicted to take 
the opposite sign. Alternatively, for an explanation related to some 
element of the AGN-fueling or life-cycle, individual quasars could 
experience episodes where both the luminosity and the E{B — V) 
increase, causing them to migrate brightward in the quasar LF for 
a period determined by the lifetime of the fueling cycle. 

While the number of quasars in our Faint Sample is modest, 
if we calculate the observed number density of reddened quasars 
down to K =18.9 (Bright) and 19.3 (Faint) we obtain surface den¬ 
sities of 0.04±0.01 deg“^ and 0.04±0.03 deg“^ respectively. The 
Poisson error associated with the Faint Sample is large but the re¬ 
sult provides additional evidence that the LF of the reddened quasar 
populations continues to increase only slowly down to fainter ab¬ 
solute magnitudes. 

lAssef et ^ l l2014h came to the same conclusion using an in¬ 
dependently selected sample of more highly obscured AGN from 
WISE. Our sample covers a larger dynamic range in terms of in¬ 
trinsic quasar luminosity as well as selecting quasars that fall at the 
lower boundary of their i?( 73-1/) distribution. In combination, the 
two investigations provide evidence for the existence of a signifi¬ 
cant population of optically-obscured, hyperluminous quasars with 
E{B — V) > 0.5, with our results probing smaller reddening s 
and fainter intrinsic luminosities than those in lAssef et alj (I 2 OI 4 I) . 
Quasars with E{B — H)>L5 are not included in our sample as 
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their observed A'-band magnitudes put t hem below the flux limit 
of our survey. At the bright-end however, lAssef et ^ ( 1201 4) have 
confirmed the existence of a population of quasars with such ex¬ 
tinctions. The space-densities shown in Fig.j^therefore represent a 
lower limit to the numbers of obscured quasars that may be missing 
from optical surveys. 


7 CONCLUSIONS 

We have presented a search for heavily dust-obscured luminous 
quasars at 2 > 2 selected from the new wide-held infrared pho¬ 
tometry provided by the UKIDSS Large Area Survey, the ESO 
VISTA Hemisphere Survey and the WISE All Sky Survey. From 48 
candidates, we used VLT-SINFONI NIR spectroscopy to conhrm 
24 new 2 > 2, heavily reddened broad-line quasars with 0.5 < 
E{B — U) < 1.5. Three of these new quasars (VHSJ2028—5740, 
ULASJ2315+0143 and ULASJ0123-f 1525) have E{B - V) > 1 
and correspond to the dustiest broad-line Type 1 quasars known. 
Combining the new quasars with objects presented in B12 & B13, 
the full sample consists of 38 spectroscopically conhrmed, heavily 
reddened quasars, 36 of which are at 2 > 2. There is no evidence 
for these near infrared selected reddened quasars having systemat¬ 
ically different Ha equivalent widths relative to their unobscured 
counterparts. We present a space density estimate for the reddened 
quasar population at 2 ~ 2 — 3 enabling detailed comparisons to 
the space-density of unobscured quasars over the same redshift and 
luminosity range. In particular, we reach the following conclusions: 

• The reddened quasars have very high bolometric lumi¬ 
nosities and black-hole masses that are comparable to the 
brightest unobscured quasars at these redshifts. The median 
bolometric luminosity of the new sample presented here is 
logio(Lboi/erg s“^)=47.1±0.4 and the median black-hole mass is 
logio(MBH/M0)=9.7±O.4, consistent with our previous samples of 
reddened quasars. There is some evidence for the presence of a tail 
of quasars out to very high bolometric luminosities of ~ 10"*® erg/s. 
Around 70 per cent of the reddened quasars have Eddington ratios 
of i/Z/Edd >0.1 confirming the presence of massive black-holes 
accreting at a high rate. Apart from the high dust extinctions along 
the line-of-sight, the properties of the reddened quasars appear to 
be similar to luminous unobscured quasars. 

• We compile WISE photometry out to 22 pm for all our red¬ 
dened quasars and using a simple power-law fit to the rest-frame 
1-4 pm SED, we show that the power law indices in the reddened 
quasars are consistent with those measured for unobscured quasars 
from the SDSS over the same redshift range and down to the same 
ZT-band magnitude limit. The mean power-law index (A Fa oc 

A/^NIR) 

at rest-frame wavelengths of l-4pm is /3nir=0.69±0.29 for 
the reddened quasars compared to 0.50±0.26 for the unobscured 
quasars. The reddened quasars have slightly steeper power-law in¬ 
dices on average due to the non-negligible extinction at rest-frame 
wavelengths of 1 — 2 pm for some of our reddest objects. 

• We use the subset of 25 22 pm-detected reddened quasars 
to calculate their rest-frame mid infrared 6 pm luminosity. We 
find an average value of logiofLe^im/erg s“^)=47.1±0.4 com¬ 
pared to logiofLe^m/erg s“^)=46.8±0.2 for equivalent SDSS 
22 pm-detected quasars. The most luminous WISE HyLIRG, 
WISEJl814-1-3412 has logio(L6^,m/erg s"i)=47.30±0.05 and we 
find seven of our sample are more luminous at 6 pm than 
WISEJ1814-I-3412. The 6 pm rest-frame luminosities are robust, 
estimated purely from the observed WISE photometry without as¬ 
sumptions regarding the quasar SED or the extinction law, corrob¬ 


orating our conclusion that these reddened quasars constitute some 
of the most luminous quasars known at these epochs. 

• We present the first estimate of the space density of reddened 
(0.5 < E{B — U) < 1.5) quasars at 2 ~ 2 — 3. We find that the 
most luminous reddened quasars, with de-reddened Mi < —29, 
are at least as numerous at these epochs as unobscured quasars. We 
also find however that the reddened quasar luminosity function has 
a much flatter slope and the number density of unobscured quasars, 
due to the rapidly rising luminosity function, greatly exceeds that 
of the reddened population by Mi ~ —27.5. The slower increase in 
the number of reddened quasars as a function of decreasing lumi¬ 
nosity suggests that as a whole, luminous reddened Type 1 quasars 
likely do not account for a significant fraction of the total energy 
liberated by the Type 1 AGN population in an ‘obscured’ phase. 

• The observed evolution of of the reddened quasar space- 
density with luminosity is also opposite to that seen in the much 
less luminous Type 2 AGN samples and to the predictions from the 
receding torus model, where the obscured fraction is predicted to 
increase with decreasing luminosity. Our population of reddened 
quasars therefore probe a very different region of the luminosity- 
extinction parameter space compared to these Type 2 AGN sam¬ 
ples, which have previously been selected through deep, small-area 
surveys in the X-ray and mid infra-red. 

• We test whether different assumptions for the reddened quasar 
luminosity function, based on the luminosity function of unob¬ 
scured quasars, are able to reproduce our observations. A model 
where the luminosity and density are evolving with redshift (LEDE 
model), and the normalisation is set by the space-density of un¬ 
obscured quasars, produces too few reddened quasars if the ex¬ 
tinction is independent of luminosity. A pure luminosity evolution 
(RLE) model with a brighter break luminosity reproduces better 
the observed counts of reddened quasars in the brightest absolute 
magnitudes in our survey, but in general results in too many red¬ 
dened quasars being predicted in our survey. Neither model to¬ 
gether with the assumption that extinction is independent of lu¬ 
minosity, is able to reproduce the observed number counts of red¬ 
dened quasars which demonstrates that either reddened quasars do 
not have the same intrinsic luminosity distribution as unobscured 
quasars, or that the extinction must depend on luminosity. A more 
elaborate model, involving an increase in the probability of finding 
high extinction quasars at high luminosities, would be consistent 
with the observations. Qualitatively, this behaviour could be related 
to some element of the AGN fuelling cycle where both the extinc¬ 
tion and luminosity increase for a period, causing the quasars to 
migrate brightward in the luminosity function. 

We conclude that our sample of high redshift (2 ~ 2.5), heav¬ 
ily reddened (0.5< E{B — V) < 2), hyperluminous (Lboi > 
10^®L©) quasars effectively bridge the gap between the optically- 
luminous, unobscured, quasars known at these redshifts and the 
much more heavily obscured AGN (E{B — U) > 5) now emerg¬ 
ing from the WISE All-Sky Survey. We demonstrate that sensitive 
near-infrared observations such as those provided by the UKIDSS 
Large Area Survey (in the northern hemisphere) and the VISTA 
Hemisphere Survey (in the southern hemisphere), are critical for 
identifying the population of quasars with intermediate extinctions 
at high redshift. The space density of such obscured quasars dom¬ 
inates the total quasar counts at very bright luminosities but they 
are unlikely to account for a significant fraction of the total en¬ 
ergy liberated by the luminous. Type 1 AGN population. Reddened 
quasars likely constitute a brief evolutionary phase in galaxy for¬ 
mation associated with high luminosities and moderate extinctions 
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Figure 9. Completeness-corrected, VIVmax estimates of the cumulative number density of reddened quasars (red solid) compared to the equivalent cumulative 
numbe r-densities of unobscured quasars derived over the s ame redshift and l uminosity range assuming the SDSS+BOSS LEDE LF model from iRoss et alJ 
iioil (blue dashed) and the SDSS+BOSS PLE LF from l^ss et al (cyan dot-dashed). Reddened quasars appear to dominate the counts at high 

intrinsic luminosities but their number counts flatten off relative to unobscured quasars as we approach more typical quasar luminosities. 


towards the quasar line-of-sight. Discriminating between different 
models for the evolution of these AGN, as they experience changes 
in their black hole growth rate, L/L^dd, degree of obscuration and 
outflow properties, will require multi-wavelength investigations of 
carefully selected sub-samples of objects and is now becoming pos¬ 
sible with sensitive NIR integral field units (e.g. SINFONI on the 
VLT) and ALMA. 
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APPENDIX A: SPECTRA & SPECTRAL LINE FITS 

We present the SINFONI JT-band spectra for all our spectroscopi¬ 
cally confirmed reddened quasars along with the best-fit SED mod¬ 
els described in Section ldTI in Fig. lAll We note that the SED mod¬ 
els include a scaling of the Ha line to match the observed equiv¬ 
alent width in each quasar. However, no effort has been made to 
match the velocity widths of the lines as these do not affect the red¬ 
dening estimates derived from this SED-fitting. The Ha line-fits 
are shown in Fig. IA2l in the velocity range ±10000 km/s. Where 
double Gaussian fits were possible as described in Section |4^ the 
fits show these double Gaussians. For the remaining objects, the fits 
correspond to a single broad Gaussian. The wavelength dependent 
noise spectrum is shown as the dashed line in the upper panel while 
the lower panel shows the fit residuals. The rest-frame full-width- 
half-maxima for both the single Gaussian and double Gaussian fits, 
are given in Table lAl] 
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VHSJ2024-5623 z=2.282 VHSJ2028-4631 z=2.464 




VHSJ2028-5740 z=2.121 VHSJ2048-4644 z=2.182 




VHSJ2100-5820 z=2.36 


VHSJ2101-5943 z=2.313 






Figure Al. VLT-SINFONI JC-band spectra of spectroscopically confirmed reddened quasars (solid line) together with our best-fit SED model (dashed line). 
Note that the equivalent width of the Ha line in the SED model has been adjusted to match the observed equivalent width of each quasar. However, all the 
SED models assume the same velocity width for the Ha line and no effort has been made to match the velocity widths as these do not affect the reddening 
estimates derived from the SED-fitting. 
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Table Al. Rest-frame FWHM from both single and double Gaussian fits to the Hq line as well as the rest-frame Hq equivalent width (EW) for each of our 
reddened quasars. The double gaussian fit is only considered if the individual components in this Gaussian have FWHM between 1000 and 10000 km/s. 


Name 

FWHMsingie/km s ^ 

FWHMdoublel/km s ^ 

FWHMdouble2 / km s ^ 

Rest-frame Ha EW / A 

VHSJ2024-5623 

7000 

1600 

8700 

-350 

VHSJ2028-4631 

4300 

- 

- 

-200 

VHSJ2028-5740 

6100 

- 

- 

-290 

VHSJ2048-4644 

4000 

- 

- 

-250 

VHSJ2100-5820 

3100 

- 

- 

-200 

VHSJ2101-5943 

5800 

3300 

9600 

-360 

VHSJ2115-5913 

3000 

- 

- 

-170 

VHSJ2130-4930 

4400 

- 

- 

-130 

VHSJ2141-4816 

3900 

1100 

4600 

-210 

VHSJ2212-4624 

5400 

3100 

7900 

-490 

VHSJ2220-5618 

3800 

1600 

5200 

-140 

VHSJ2227-5203 

7300 

- 

- 

-190 

VHSJ2235-5750 

4900 

2100 

9000 

-260 

VHSJ2256-4800 

3900 

2600 

8800 

-200 

VHSJ2257-4700 

4000 

2100 

5200 

-420 

VHSJ2306-5447 

7900 

1500 

8600 

-500 

VHSJ2332-5240 

5900 

- 

- 

-300 

VHSJ1556-0835 

1900 

- 

- 

-230 

VHSJ2143-0643 

3300 

2000 

8300 

-250 

VHSJ2144-0523 

3200 

2300 

9600 

-250 

VHSJ2109-0026 

3500 

1900 

6900 

-320 

VHSJ2355-0011 

6000 

2400 

9000 

-330 

ULASJ0123-H1525 

3500 

1000 

4100 

-380 

ULASJ2315-H0143 

4700 

3000 

7000 

-270 
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Figure A1 - continued 
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Figure Al - continued 
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VHSJ2143-0643 z=2.383 



Observed A. (|i m) 
VHSJ2109-0026 z=2.344 



Observed X (}i m) 
ULASJ0123+1525 z=2.629 



VHSJ1556-0835 z=2.188 



VHSJ2144-0523 z=2.152 



Observed X (fi m) 


VHSJ2355-0011 z=2.531 



Observed X (}i m) 


ULASJ2315+0143 z=2.56 



Observed X (}i m) 


© 0000 RAS, MNRAS 000, 000-000 

















































































Heavily Reddened Type 1 Quasars at z > 2 21 


VHSJ2024-5623 z=2.282 E(B-V)=0.6 
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Figure A2. Hq line profiles and Gaussian line fits. Bottom panel shows the residuals and the dashed line is the wavelength dependent noise spectrum. 
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Figure A2 - continued 


VHSJ2141-4816 z=2.655 E(B-V)=0.8 



VHSJ2212-4624 z=2.141 E(B-V)=0.8 
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Table Bl. K-corrections in the i-band and i^-band derived using our quasar 
SED model from Section|4J] 


Redshift 

ki 

kK 

2.00 

-1.084 

0.506 

2.10 

-1.078 

0.287 

2.20 

-1.069 

0.151 

2.30 

-1.061 

0.118 

2.40 

-1.064 

0.091 

2.50 

-1.071 

0.092 

2.60 

-1.099 

0.131 

2.70 

-1.153 

0.333 

2.80 

-1.181 

0.345 

2.90 

-1.195 

0.321 

3.00 

-1.196 

0.292 


APPENDIX B: SED MODEL AND K-CORRECTIONS 

Our computations of the quasar bolometric luminosity, absolute 
magnitude and resulting space density relative to unobscured 
quasars, implicitly assume the quasar SED model described in Sec¬ 
tion |4T] As this SED model is different from several others used 
in the literature, we provide in Table HD the i-band and if-band 
k-corrections for this model over the redshift range of our survey. 
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Figure A2 - continued 


VHSJ2332-5240 z=2.45 E(B-V)=0.6 



VHSJ1556-0835 z=2.188 E(B-V)=0.7 
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